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Reynolds-number criteria are developed for acceptable variations in Space Shuttle Or- 
biter entry trajectories for use in computational aeroheating analyses. The criteria deter- 
mine if an existing computational fluid dynamics solution for a particular trajectory can 
be extrapolated to a different trajectory. The criteria development considers twelve types 
of computational aeroheating data, such as boundary layer thickness. For each type of da- 
tum, the allowable uncertainty contribution due to trajectory variation has been set by the 
Entry Aeroheating Subsystem team. Then Reynolds-number relations between trajectory 
variation and output uncertainty are determined. From these relations the criteria are es- 
tablished for the maximum allowable trajectory variations. The most restrictive criterion 
allows a 25% variation in Reynolds number at constant Mach number between trajectories. 

I. Nomenclature 

The subject matter in the present report considers relative changes in fluid dynamic properties. Therefore 
dimensional units for the symbols are not pertinent to the results. 

Symbols 

Cp Pressure coefficient 

D Cavity depth 

H Total enthalpy 

M Mach number 

P Pressure 

q Heat transfer rate 

Re Reynolds number 

T Temperature 

V Velocity magnitude 

6 Boundary layer thickness 

S* Displacement thickness 

H Viscosity 

p Density 

9 Momentum thickness 


Subscripts 

e Boundary layer edge 

k Protuberance height 


* Aerothermodynamics Branch, senior AIAA member. 

* Senior Research Scientist, AIAA Associate Fellow. 
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w 


Wall 

Free stream 


oo 


Operators 

A Difference: A a: = X 2 — aq 


Acronyms 


BF 

BLT 

CFD 

DPLR 

LAURA 

OML 

RCC 


Bump factor 

Boundary layer transition 
Computational fluid dynamics 
Data-parallel line relaxation CFD code 

Langley aerothermodynamic upwind relaxation algorithm CFD code 
Outer mold line 
Reinforced carbon-carbon 


II. Introduction 

The Space Shuttle Orbiter Entry Aeroheating Subsystem team has pre-computed computational fluid 
dynamics (CFD) solutions of Orbiter flow fields for nominal entry trajectories. 1 During missions, these 
solutions may be used to predict aerothermodynamic properties in the course of assessing Orbiter damage 
and repairs. 2 But, the mission-specific entry trajectory will not match the nominal trajectory, because of 
differences such as the Orbiter weight or atmospheric conditions. These differences in trajectories can be 
expressed as differences in Reynolds numbers at given Mach numbers. This report develops Reynolds- number 
criteria to characterize how close in trajectory space a particular trajectory must be to a nominal trajectory 
so that the nominal CFD solutions are sufficient for damage assessment purposes. The most restrictive 
criterion is then used to define the allowable trajectory variations for reuse of the nominal CFD solutions 
during a mission. 

The criteria are developed under several assumptions: a spot-check at one Mach number applies for all 
Mach numbers (the Mach numbers of most interest are greater than 15); uncertainties are symmetric and 
can be expressed as ± quantities; uncertainty trends extrapolate linearly in Reynolds number; and, angle of 
attack variations are negligible (the angles of attack are within 40 ± 2°). 

The development of the criteria begins by identifying twelve types of aerothermodynamic data, such as 
surface temperature and boundary layer thickness, that are available in the CFD solution set and are relevant 
to Orbiter damage or repair assessments. For each datum type, accuracy requirements at approximately 
1:20 odds have been set by the Entry Aeroheating Subsystem team. Reynolds number scalings are then 
developed to relate trajectory variations to variations in the CFD results. Finally, the allowable trajectory 
Reynolds number variations are determined by inverting the scaling relations. 

The Reynolds number criteria are summarized as a tabulation in section IV-M. 

An expanded treatment of the present material, along with similar analyses for an additional five more 
types of aerothermodynamic data, is available in a NASA Technical Memorandum. 3 

III. Nominal Orbiter CFD solutions 

Aeroheating CFD solutions of entry flow fields about the Orbiter vehicle 3, exist for several points along a 
few nominal trajectories. The CFD codes, the solution processes, the verification techniques, the validation 
data, and the solutions themselves have been accepted for use during Space Shuttle missions by the Orbiter 
Configuration Control Board. 1 Previous Orbiter solutions of a similar nature and application, using the 
same CFD codes, have been reported 4 in the context of the Columbia accident investigation. 

These solutions typically require about a week each to obtain, although with exclusive access to a super- 
computer in a round-the-clock staffing mode solutions have been obtained in a day or two. It is desired to 

a The CFD Orbiter model used by the Entry Aeroheating Subsystem team is truncated in the vicinity of the body-flap hinge 
line, 1280 inches from the nose. 
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Priority 
1 
2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

Table 1. CFD aeroheating products, ordered by perceived importance. 

avoid this cost of creating new entry Orbiter solutions; thus the present report’s topic of how much these 
existing solutions can be reused for different trajectories. 

A. CFD codes 

Two CFD codes are normally employed by the Entry Aeroheating Subsystem team, and the present analyses 
utilize the solutions from both codes interchangeably. The two codes provide a measure of verification when 
used redundantly, and pragmatically allows the team members to choose their most proficient analysis tool. 
The two codes are the Data Parallel Line-Relaxation (DPLR) software 5 and the Langley Aerothermodynamic 
Upwind Relaxation Algorithm (LAURA ). 6, 7 

Both DPLR and LAURA are finite- volume upwind-flux schemes for solving the Navier-Stokes 8,9 equations 
with non-equilibrium chemistry. All the CFD results considered here have laminar boundary layers, a 
five-species air model, and assume a radiative-equilibrium surface temperature using reaction-cured glass 
properties. 

B. Solution repository 

There are currently 20 full-vehicle Orbiter solutions covering parts of three nominal trajectories. The so- 
lutions are stored on a data repository at the NASA Advanced Supercomputing facility. The CFD data 
for the present report was all obtained from this solution repository. The solutions span the Mach number 
range 6-25. The trajectories correspond to STS-107 ( Columbia’s last flight), STS-115 (a typical entry from 
the International Space Station), and ISSHVFW (a severe entry from the International Space Station). 

IV. Trajectory variation criteria 

Twelve CFD aeroheating products are considered. The products are listed in table 1 as prioritized by the 
Entry Aeroheating Subsystem team, based on the perceived importance for mission-support assessments. 
The products are considered by this report only on the wind side of the Orbiter. 

For each product, the following sections define the product, explain its application for mission support, 
and present an accuracy requirement for the CFD prediction from the Entry Aeroheating Subsystem team. 
The Reynolds number is used to parameterize the degree a trajectory is off-nominal, and criteria are developed 
that relate trajectory Reynolds number variation to the allowable error in the aeroheating product due to 
trajectory variation from nominal. The criteria limits are set in the context of approximately 1:20 odds. 
These criteria are developed subject to several assumptions. Usually, spot-checks are only performed at one 
Mach number, with the criteria assumed to apply for all Mach numbers; the Mach numbers of most interest 
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Parameter 

Idealized cavity heating bump factors 
Smooth OML surface temperatures 
Smooth OML heating rates 
Boundary layer thickness 
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are 15-25. Variations are most often only considered in one direction, and the rates of change are assumed 
to be the same in both the positive and negative sides. The product variations are assumed to extrapolate 
linearly with Reynolds number. The effects of angle of attack variations are assumed to be negligible; the 
angles of attack are within 40 ± 2°. 

The trajectory criteria are summarized in table 3 in section IV-M. 


A. Cavity heating bump factors 

The cavity heating bump factors are surface heat transfer rates that have been normalized to a reference 
nominal heating rate on the smooth OML. The cavities are simplified models of actual damage geometries, 
such as hexahedrons, and solutions are obtained as local perturbations within the global Orbiter flow field. 
In the terminology of the Entry Aeroheating Subsystem team, these geometry models are called idealized 
cavities. Examples of the cavity modeling and analysis process have been published by Pulsonetti and 
Wood. 10 CFD cavity heating bump factors can be used during missions to assess and supplement the 
predictions from the Cavity Heating Tool, 11 which is the standard cavity heating engineering tool used by 
the Entry Aeroheating Subsystem team. Of primary interest is the predicted bump factors on the floor of a 
cavity. The Entry Aeroheating Subsystem team has set the allowable error in bump factors due to trajectory 
variations at ±25%. 

The trajectory criteria for cavity bump factors is developed using the Cavity Heating Tool correlations, 
which are based upon cavity wind tunnel data. In particular, two of the Cavity Heating Tool correlations 
show a sensitivity to free stream conditions: the fore-85% floor of Everhart (short) cavities and the aft-10% 
floor of long cavities. See Hyatt et al. 11 for a full description of the Cavity Heating Tool nomenclature. 

The Everhart cavity fore-85% floor bump factor varies with boundary layer thickness as 


BF = 


5 
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The difference in bump factor for the same cavity at two different boundary layer thicknesses can be expressed 
in absolute and relative terms as 


ABF = BF 2 - BFi 
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Similarly, the long cavity aft-10% floor bump factor varies with boundary layer thickness as 
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The absolute and relative changes in bump factor as functions of boundary layer thickness are 
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Variations in boundary layer thickness can be related to variations in Reynolds number, because the 
boundary layer thickness varies approximately linearly with the reciprocal of the square root of the Reynolds 
number, as demonstrated in section D. This relationship is expressed as 
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Now, the change in bump factors can be expressed in terms of Reynolds number variations, and so can 
the inverse relations for Reynolds number change be expressed in terms of bump factor variations. For the 
Everhart-cavity fore-85% floor bump factor 


ABF = - 
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For the long-cavity aft- 10% floor bump factor 
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Because equations 11, 13, 15, and 17 are monotonic in 5i, a positive-definite quantity, the extreme value 
theorem states that the maximum variation will occur at an endpoint of the <5 ranges listed in equations 1 
and 5. 

By equation 1, BF(D/S = 0.1) = 1. Using either equation 11 or 13, it can be determined that a 300% 
increase in Reynolds number is needed to reduce the bump factor by 25% to 0.75. 

By equation 1, BF(D/S = 0.3) = 0.5. Using either equation 11 or 13, it can be determined that a 31% 
decrease in Reynolds number is needed to increase the bump factor by 25% to 0.625. 

By equation 5, BF(D/S = 0.3) = 0.8. Using either equation 15 or 17, it can be determined that a 31% 
increase in Reynolds number is needed to increase the bump factor by 25% to 1.0. 

By equation 5, BF(D/S = 1.2) = 5. Using either equation 15 or 17, it can be determined that a 40% 
decrease in Reynolds number is needed to decrease the bump factor by 25% to 3.75. 

Coincidentally, the limiting allowable variation in Reynolds number to produce a 25% variation in bump 
factor is the same for both the Everhart and long cavities — 31%. Thus the cavity heating bump factor 
criteria for allowable trajectory variation is a ±31%i variation in Reynolds number. 


B. Surface temperature 

The CFD surface temperatures are radiative equilibrium values assuming reaction-cured glass properties for 
the entire Orbiter. Of primary interest is the wind-side tile acreage of the vehicle. The CFD surface temper- 
atures may be used during a mission to assess the predicted values from the standard entry aeroheating tools. 
The Entry Aeroheating Subsystem team has set the allowable change in surface temperature predictions due 
to trajectory variations at ±3%. 

The variation in surface temperature with free stream conditions has been approximately scaled by 
O^O 1 / 8 for an Orbiter boundary layer probe tool. 12 But the present analysis is considering variations in 
Reynolds number at a fixed Mach number. For the Orbiter entry trajectories, because the temperature does 
not vary much in the mesosphere, the possible velocity variation at a fixed Mach number is small. Thus, 
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the Reynolds number variations are caused primarily due to density variations at different altitudes, and the 
surface temperature changes can be related to Reynolds number variations scaled to the 1/8 power. 

The Mach-18 STS-115 trajectory point has a 21.1% lower Reynolds number than the Mach-18 STS-107 
trajectory point. Multiplying the Mach-18 STS-107 nominal surface temperatures by the Reynolds number 
ratio to the 1/8 power, 0.789^ 1//8 ^ = 0.971, produces an estimate of the STS-115 surface temperatures. The 
error in this estimate relative to the nominal STS-115 solution, for the wind side of the Orbiter, is shown in 
figure 1. The errors are within ±2% over at least 95% of the surface. 



Figure 1. Extrapolation error in surface temperature for Mach-18 check case. 


Using a ±2% bound 
surface temperatures is 


on the error, the Reynolds number criterion for allowable trajectory variation of 


/ 3% allowable \ 
\ 2% observed / 


31.7% 


±32% 


(18) 


C. Smooth OML heating 

The CFD surface heat transfer rates are based upon reaction-cured glass properties for the entire Orbiter. 
Of primary interest is the wind-side tile acreage of the vehicle. The CFD smooth-OML heating rates may be 
used during a mission to assess the predicted values from the standard entry aeroheating tools. The Entry 
Aeroheating Subsystem team has set the allowable change in heat transfer rate predictions due to trajectory 
variations at ±10%. 

In section B, the surface temperatures are shown to scale with For radiative equilibrium the 

relation between surface temperature and heating rates is q oc T 4 . The assumed variation in heating due to 
Reynolds number changes is then q oc \J Re. 

The Mach-18 STS-115 trajectory point has a 21.1% lower Reynolds number than the Mach-18 STS- 
107 trajectory point. Multiplying the Mach-18 STS-107 nominal surface heating rates by v / R789 = 0.888 
produces an estimate of the STS-115 heating rates. The error in this estimate relative to the nominal STS- 
115 solution, for the wind side of the Orbiter, is shown in figure 2. The errors are within ±7% over about 
95% of the surface. 

Using a ±7% bound on the error, the Reynolds number criterion for allowable trajectory variation of 
smooth-OML surface heating rates is 


y ) 21.1% = 30.1% ~ ±30% 


(19) 
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Figure 2. Extrapolation error in heating rates for Mach-18 check case. 


D. Boundary layer thickness 

The CFD boundary layer thickness has been defined by the Entry Aeroheating Subsystem team for the 
Orbiter flow fields as a variation on the traditional 99.5% H/Hoo definition. See references 1, 12,13 for details 
of the boundary layer thickness definition. The boundary layer thickness may be used as a correlating 
parameter when creating engineering models of Orbiter damage during a mission. The Entry Aeroheating 
Subsystem team has set the allowable change in boundary layer thickness predictions due to trajectory 
variations at ±11%. 

Classic incompressible boundary layer theory 14 scales boundary layer thickness inversely with \J Re. A 
1/y/Re scaling was used to extrapolate boundary layer thicknesses at three locations on the Orbiter — at 
the nose landing gear door, at the external tank door, and at the RCC-9 carrier panel-from the STS-107 
trajectory to the STS-115 trajectory for six trajectory points spanning Mach 17.9-24.9. The scaling factors 
are listed in table 2. 


Moo 

Re 115 Re 107 

/ Re 107 

Re 107 

V Re ii 5 

24.9 

19.6% 

0.91 

24.2 

6.2 

0.97 

22.9 

- 9.4 

1.05 

20.3 

- 7.1 

1.04 

19.4 

-16.7 

1.10 

17.9 

-21.1 

1.13 


Table 2. Boundary layer thickness scaling factors from STS-107 to STS-115 trajectories. 


Figure 3 plots the absolute value of the extrapolation error, relative to CFD solutions for the STS-115 
free-stream conditions, for each of the eighteen check points versus the absolute value of the percentage 
variation in Reynolds number between the STS-107 and STS-115 trajectories. No simple correlation is 
evident, but the errors are bounded by about 5.5%, with a maximum 21% variation in Reynolds number. 

The Reynolds number criterion for allowable trajectory variation of boundary layer thicknesses is 

^ x 21.1% = 42.2% ~ ±40% (20) 
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Figure 3. Extrapolation errors for boundary layer thickness test points. 


E. Protuberance boundary layer trip parameter 

As part of the primary correlation for predicting boundary layer transition induced by a protuberance, the 
Entry Aeroheating Subsystem team uses the parameter 


Ree 1 

Me ' 6 ' \H W J 


(21) 


All of the variables in parameter (21) are obtained from CFD solutions of the smooth-OML Orbiter. The 
scaling of this parameter leads to the most restrictive Reynolds number trajectory criterion. The Entry 
Aeroheating Subsystem team has set the allowable error contribution due to trajectory variations at ±10%. 

The boundary layer scalings from the Boundary Layer Transition Tool documentation 12 are used to 
obtain a scaling for parameter (21): M e is not scaled with Re, Reg is scaled with \J~Re, and <5 is scaled with 
l/VKe. (H e /H w ) 03 is shown in reference 12 to scale with Re/10, and is neglected — meaning it is unsealed 
in the present analysis — because the variation is an order of magnitude less than the variation of Reg/5. 
Thus the parameter (21) is scaled with Re. 

The linear scaling with Reynolds number is tested at Mach 19.4, by extrapolating the STS-107 parameter 
to the STS-115 conditions, a 16.7% reduction in Reynolds number. Absolute values of the percentage 
differences between the computed and the extrapolated STS-115 parameter-(21) values that exceed 6.5% 
are shown as red in figure 4. Approximately 95% of the windside surface has extrapolation errors less than 
6.5%. 

The Reynolds number criterion for the allowable trajectory variation of the protuberance boundary layer 
trip parameter is 

^ x 16.7% = 25.7% ~ ±25% (22) 
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Figure 4. Absolute values of percentage extrapolation errors for protuberance boundary layer transition 
parameter (21). Red is more than 6.5%. Blue is less than 6.5%. 



Figure 5. Absolute values of percentage extrapolation errors for Re^ boundary layer transition parameter (23). 
Red is more than 4.5%. Blue is less than 4.5%. 


F. Rek boundary layer transition parameter 

As part of the secondary correlation for predicting boundary layer transition induced by a protuberance, the 
Entry Aeroheating Subsystem team uses the parameter 

ite°- 6 (W^) (23) 

V Mfc/ 

All of the variables in parameter (23) are obtained from CFD solutions of the smooth-OML Orbiter. The 
Entry Aeroheating Subsystem team has set the allowable error contribution due to trajectory variations 
at ±10%. 

Assuming k ~ <5, then both Rek and Reg can be scaled with %/ Re. Viscosity scales approximately with 
the square root of temperature, and so ( /i e /fik) 0A scales like (T e /Tk) 0 ' 2 , and is neglected. Substitution into 
the expression (23) leads to the entire parameter being scaled with \J~Re. 

The \[Re scaling is tested at Mach 19.4, by extrapolating the STS-107 parameter (23) to the STS-115 
conditions, a 16.7% reduction in Reynolds number. Absolute values of the percentage differences between 
the computed and the extrapolated STS-115 parameter- (23) values that exceed 4.5% are shown as red in 
figure 5. Approximately 95% of the windside surface has extrapolation errors less than 4.5%. 

The Reynolds number criterion for the allowable trajectory variation of the Rek boundary layer trip 
parameter is 

^ x 16.7% = 37.1% ~ ±37% (24) 
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Figure 6. Absolute values of percentage extrapolation errors for cavity-volume boundary layer transition 
parameter (25). Red is more than 4%. Blue is less than 4%. 


G. Cavity-volume boundary layer transition parameter 


As part of the volume-based correlation for predicting boundary layer transition induced by a cavity, the 
Entry Aeroheating Subsystem team uses the parameter 


Reg 1 
M e ' S 



(25) 


All of the variables in parameter (25) are obtained from CFD solutions of the smooth-OML Orbiter. The 
Entry Aeroheating Subsystem team has set the allowable error contribution due to trajectory variations 
at ±10%. 

Reg is scaled with yRe, S is scaled inversely with yRe, and M e is unsealed. The temperature ratio 
scaling is neglected, because T e /T w is shown to scale with Re 018 in reference, 12 and when T e /T w is raised to 
the 0.51 power the scaling becomes less than Re 0 ' 1 . The entire parameter (25) is thus scaled linearly with 
Reynolds number. 

The linear scaling with Reynolds number is tested at Mach 19.4, by extrapolating the STS-107 parameter 
to the STS-115 conditions, a 16.7% reduction in Reynolds number. Absolute values of the percentage 
differences between the computed and the extrapolated STS-115 parameter-(25) values that exceed 4% are 
shown as red in figure 6. Approximately 95% of the windside surface has extrapolation errors less than 4%. 

The Reynolds number criterion for the allowable trajectory variation of the cavity-volume boundary layer 
trip parameter is 

^ x 16.7% = 41.8% ~ ±40% (26) 



H. Cavity-area boundary layer transition parameter 


As part of the area-based correlation for predicting boundary layer transition induced by a cavity, the Entry 
Aeroheating Subsystem team uses the parameter 


Reg 1 

~M~e ' 6 



(27) 


All of the variables in parameter (27) are obtained from CFD solutions of the smooth-OML Orbiter. The 
Entry Aeroheating Subsystem team has set the allowable error contribution due to trajectory variations 
at ±10%. 

Reg is scaled with yRe, 5 is scaled inversely with yRe, and M e is unsealed. The temperature ratio 
scaling is neglected because T e /T w is shown to scale with Re 018 in reference, 12 and when T e /T w is raised to 
the 0.67 power the scaling becomes approximately Re 0 ' 1 . The entire parameter (27) is thus scaled linearly 
with Reynolds number. 

The linear scaling with Reynolds number is tested at Mach 19.4, by extrapolating the STS-107 parameter 
to the STS-115 conditions, a 16.7% reduction in Reynolds number. Absolute values of the percentage 
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Figure 7. Absolute values of percentage extrapolation errors for cavity-area boundary layer transition param- 
eter (27). Red is more than 4%. Blue is less than 4 %. 


Figure 8. Absolute values of percentage errors for post normal-shock stagnation pressure coefficients. Black 
is more than 3 %. 

differences between the computed and the extrapolated STS-115 parameter-(27) values that exceed 4% are 
shown as red in figure 7. Approximately 95% of the windside surface has extrapolation errors less than 4%. 


I. Boundary layer C p , shock relations 

One of the techniques used by the Entry Aeroheating Subsystem team to estimate protuberance bending loads 
utilizes boundary-layer profile pressure coefficients that assume normal-shock stagnation pressures. These 
pressure coefficients are obtained from the smooth-OML CFD baseline solutions. The Entry Aeroheating 
Subsystem team has set the allowable error contribution due to trajectory variations at ±12%. 

To test the boundary layer pressure coefficient scaling, post normal-shock pressure coefficients are com- 
puted at the boundary layer edge for Mach-19.4 trajectory points on both the STS-107 and STS-115 trajec- 
tories. These two trajectory points differ by 16.7% in Reynolds number. 

The absolute values of the percentage differences between the wind-side STS-115 and STS-107 pressure 
coefficients are shown in figure 8. Black regions have differences greater than 3%. Approximately 95% of 
the wind-side surface has differences less than 3%. 

The Reynolds number criterion for the allowable trajectory variation of the normal-shock boundary layer 
pressure coefficients is 





( 29 ) 


11 of 15 


American Institute of Aeronautics and Astronautics 



Figure 9. Absolute values of percentage errors for total-pressure coefficients. Black is more than 3%. 


J. Boundary layer C p , stream-tube 

One of the techniques used by the Entry Aeroheating Subsystem team to estimate protuberance bending 
loads utilizes boundary-layer profile pressure coefficients that assume stream-tube theory total pressures, 
P ± \pV 2 . These pressure coefficients are obtained from the smooth-OML CFD baseline solutions. The 
Entry Aeroheating Subsystem team has set the allowable error contribution due to trajectory variations 
at ±12%. 

To test the boundary layer pressure coefficient scaling, total-pressure coefficients are computed at the 
boundary layer edge for Mach-19.4 trajectory points on both the STS-107 and STS-115 trajectories. These 
two trajectory points differ by 16.7% in Reynolds number. 

The absolute values of the percentage differences between the wind-side STS-115 and STS-107 pressure 
coefficients are shown in figure 9. Black regions have differences greater than 3%. Approximately 95% of 
the wind-side surface has differences less than 3%. 

The Reynolds number criterion for the allowable trajectory variation of the stream-tube boundary layer 
pressure coefficients is 

^ x 16.7% = 66.8% ~ ±67% (30) 



K. Displacement thickness 

The Entry Aeroheating Subsystem team has developed a tool, named Endigestion, for estimating breach 
flows. During operation, the Endigestion tool integrates mass flux through the boundary layer. The flow 
properties used by Endigestion are obtained from the smooth-OML CFD baseline solutions. The Entry 
Aeroheating Subsystem team has set the allowable error contribution due to trajectory variations at ±50%. 

The displacement thickness is used to estimate the trajectory error scaling of the Endigestion mass flux. 
The displacement thickness scales with 1/y/Re. 

The STS-107 S* at M <*, = 19.4 are scaled to the STS-115 trajectory, a 16.7% Re difference. The absolute 
values of the percentage scaling errors are shown in figure 10, where the red indicates errors greater than 
±5%. Approximately 95% of the wind-side surface has scaling errors less than 5%. 

The Reynolds number criterion for the allowable trajectory variation of the displacement thickness is 

x 16.7% = 167% ~ ±165% (31) 


L. Total enthalpy ratio 

The Entry Aeroheating Subsystem team has developed a tool, named Endigestion, for estimating breach 
flows. The Endigestion tool estimates the fraction of free stream total enthalpy that is ingested into the 
breach. The flow properties used by Endigestion are obtained from the smooth-OML CFD baseline solutions. 
The Entry Aeroheating Subsystem team has set the allowable error contribution due to trajectory variations 
at ±50%. 
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Figure 10. Absolute values of percentage errors for 5* scaling. Red is more than 5%. Blue is less than 5%. 



Figure 11. Absolute values of percentage errors from scaling the percentage changes in H w /H e by -0.4 times 
the percentage change in Reynolds number. Red is more than 3%. Blue is less than 3%. 


The scaling of the total enthalpy ratio H w /H e is used here to estimate the trajectory scaling error for the 
Endigestion total enthalpy ratio. As part of the BLT (version 2) documentation, 12 the percentage change in 
the ratio H w /H e is shown to scale as —0.4 times the percentage change in Reynolds number for > 20. 
In reference, 12 this scaling has been normalized for a 10% change in Reynolds number. Figure 11 re-plots 
the scaled data from reference, 12 this time showing a threshold of ±3% error in the scaling, marked red in 
the figure. Over 95% of the wind side surface is seen to have a scaling error less than 3%. 

The Reynolds number criterion for the allowable trajectory variation of the total enthalpy ratio is 

x 10% = 167% ~ ±165% (32) 


M. Tabulation 

Table 3 lists the allowable error contributions due to trajectory variations for each of the twelve CFD 
products discussed in the present report. The last column of table 3 lists the resulting allowable trajectory 
variation, expressed in terms of trajectory Reynolds number. The most restrictive criterion, for the primary 
protuberance-induced boundary layer transition parameter, allows a ±25% variation in trajectory Reynolds 
number. 


V. Summary of results 

Criteria were established for reusing the computational fluid dynamics (CFD) solutions of the Space 
Shuttle Orbiter at specific entry trajectory points for new entry trajectories. These criteria were developed 
for twelve particular types of aeroheating data that the Entry Aeroheating Subsystem team might need during 
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Priority Parameter 


allowable trajectory trajectory Re 

error contribution, % criterion, % 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


Idealized cavity heating bump factors 
Smooth OML surface temperatures 
Smooth OML heating rates 
Boundary layer thickness 

Ree l ( g e V 
M e ' 5 ' \H W J 


0.3 


Re™ ( Re e 

*■(£) 

L(ft) 


Re». 

M. 


0.4 
He 
Hk 
0.51 


Re t 

~W. 


0.67 


Boundary layer C p , shock relations 
Boundary layer C p , stream-tube relations 
Displacement thickness 

H,„ 

He 


25 

31 

3 

32 

10 

30 

11 

40 

10 

25 

10 

37 

10 

40 

10 

40 

12 

67 

12 

67 

50 

165 

50 

165 


Table 3. Summary of allowable error contributions due to trajectory variations for CFD products and the 
resulting trajectory variation limits expressed as Reynolds number thresholds. 


a Space Shuttle mission to assess thermal protection damage to the Orbiter. The criteria are expressed in 
terms of Reynolds number variations for fixed Mach numbers. 

The development of the trajectory variation criteria began by defining the allowable error contribution 
due to trajectory variations as set by the Entry Aeroheating Subsystem team for each of the twelve CFD 
aeroheating datum types. Scaling relations were applied to cast the datum uncertainty in terms of the 
trajectory Reynolds number variations. Spot checks of the scalings were performed against CFD solutions 
for different trajectories, and the allowable trajectory variation criteria were obtained to approximately 
95%-confidence levels. 

The criteria were developed under several assumptions: a spot-check at one Mach number applied for all 
Mach numbers (the Mach numbers of most interest were greater than 15); uncertainties were the same for 
the positive and negative sides; uncertainty trends extrapolated linearly in Reynolds number; and, angle of 
attack variations were negligible (the angles of attack were within 40 ± 2°). 

The most restrictive trajectory variation criteria, for a protuberance-induced boundary layer transition 
correlation parameter, allows for a 25% variation in Reynolds number. The criteria for all the eight highest- 
priority datum types allow 40% or less variations in Reynolds number. In light of the noted development 
assumptions, for future Space Shuttle missions the Entry Aeroheating Subsystem team has decided to use 
a criteria of 15-20% variation in Reynolds number, when assessing the need for trajectory-specific CFD 
solutions. 
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